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[i] A unique microphysical structure of rainfall is observed 
by the surface laser optical Particle Size and Velocity 
(Parsivel) disdrometers on 25 April 2011 during Midlati- 
tude Continental Convective Clouds Experiment (MC3E). 
According to the systematic differences in rainfall rate and 
bulk effective droplet radius, the sampling data can be 
divided into two groups; the rainfall mostly from the deep 
convective clouds has relatively high rainfall rate and large 
bulk effective droplet radius, whereas the reverse is true for 
the rainfall from the shallow wrm clouds. The Weather 
Research and Forecasting model coupled with spectral bin 
microphysics (WRF-SBM) successfully reproduces the two 
distinct modes in the observed rainfall microphysical struc- 
ture. The results show that the up-to-date model can dem- 
onstrate how the cloud physics and the weather condition 
on the day are involved in forming the unique rainfall 
characteristic. Citation: Iguchi, T., T. Matsui, A. Tokay, P. 
Kollias, and W.-K. Tao (2012), Two distinct modes in one-day 
rainfall event during MC3E field campaign: Analyses of disdrom- 
eter observations and WRF-SBM simulation, Geophys. Res. Lett., 
39, L24805, doi:10.1029/2012GL053329. 

1. Introduction 

[2] Precipitation is a key component of earth’s energy and 
water cycle. It is an integral product of hydrometeor size 
distribution where the size of the hydrometeors ranges from 
a small drizzle drop (>0.1 mm in diameter) to a hailstone 
(>10 cm in diameter). Several processes are responsible for 
determining the size and phase of individual hydrometeors at 
the surface including but not limited to aerosol loading, 
liquid water availability, structures of vertical and horizontal 
wind velocities, presence or absence of mixed-phase pro- 
cesses, and subcloud layer evaporation. 

[ 3 ] Microphysical structure of precipitation at the ground 
is of interest in developing cloud microphysical models, in 
terms of improving forecast accuracy of surface precipitation 
rate as well as validating cloud microphysics by regarding 
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precipitation as a product. Bulk cloud microphysics is com- 
monly used in the present atmospheric models; however, 
there is little room for change in the precipitation micro- 
physics because certain mathematical functions are used as 
built-in parameterization to represent particle size distribu- 
tion (PSD). In contrast, spectral bin microphysics (SBM) 
allows prediction of hydrometeor PSDs and subsequently 
explicit representation of precipitation microphysical struc- 
ture. Numerical weather prediction (NWP) simulations using 
SBM have been utilized mostly in the research field of cloud- 
aerosol interactions in previous studies [e.g., Lynn et al., 
2005]. Analysis and verification of the ground precipitation 
microphysics can be a new application of the simulation to 
show possibilities of the up-to-date models. 

[4] This paper reports on the analysis of distinct rainfall 
microphysical structures associated with deep convective 
system and shallow boundary layer clouds in the rainfall 
event on 25 April 2011. The event was observed during the 
Midlatitude Continental Convective Clouds Experiment 
(MC3E) at the Atmospheric Radiation Measurement (ARM) 
Southern Great Plains (SGP) site in central Oklahoma. The 
disdrometer-based observations for raindrop size distribution 
(DSD) are analyzed to determine the characteristic of the 
rainfall, which is compared with a simulation result of the 
Weather Research and Forecasting model coupled with 
spectral bin microphysics (WRF-SBM). 

2. Description of Model Simulation and Overall 
Condition on 25 April 2011 at SGP Site 

[5] The WRF-SBM has been developed for supporting the 
National Aeronautics and Space Administration (NASA) 
Global Precipitation Measurement (GPM) mission; the MC3E 
field campaign is a part of GPM Ground validation (GV) 
program. The SBM part is the same as the microphysical 
scheme of the Hebrew University Cloud Model (HUCM) 
[. Khain et al., 2011]. The scheme has been implemented into 
the Advanced Research WRF (ARW) version 3.1.1 [Iguchi 
et al., 2012], Hydrometeors are categorized into one-water 
and six-ice classes, i.e., water droplets, ice crystals (plate, 
column, dendrite), snow aggregates, graupel, and hail; the 
discrete PSDs of the hydrometeor classes are represented on 
a grid containing 43 doubling mass bins. In order to simu- 
late the entire rainfall event, the simulation domain covers 
454 x 388 horizontal grids domain with 3 km horizontal 
grid spacing and inner-telescoped 642 x 498 horizontal 
grids domain with 1 km grid spacing. The inner domain 
includes the location of the ARM SGP site. The WRF-SBM 
was employed in 24-hours real-time forecasting on 25 April 
2011 using online 2-way grid nesting configuration for 
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a) ARM SGP KAZR vertical reflectivity profile, 25 April 2011 
insect-filtered reflectivity (dBZ) and cloud base estimated by Celiometer 
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c) ARM SGP central facility 
Balloon-borne sounding system 
0823 UTC 25 April 2011 



d) ARM SGP central facility 
Balloon-borne sounding system 
2039 UTC 25 April 2011 



Figure 1. Time series of (a) the Ka-band ARM Zenith Radar (KAZR) vertical reflectivity profile and (b) the rainfall rate 


measured by the laser optical Parsivel disdrometer over ARM 
sounding at Cl site at (c) 0823 UTC and (d) 2039 UTC. 


the two domains (detailed information can be found in 
auxiliary material). 1 

[6] The Ka-band ARM Zenith Radar (KAZR) observa- 
tions [e.g., Clothiaux et al., 2000; Kollias et ah, 2007] at the 
SGP central facility site show that a series of deep convective 
clouds with the cloud top heights up to 10 km passed over the 
site from 09 UTC to 1 1 UTC on 25 April 201 1 (Figure la). 
High rainfall rate up to approximately 40 mm hr -1 was 
observed by the collocated laser optical Particle Size and 
Velocity (Parsivel) disdrometer during this period (Figure lb). 
On the other hand, shallow boundary layer clouds less than 
the height of 2 km continuously hanged over the site, and 
weak rainfall approximately 1 mm hr - 1 was observed inter- 
mittently from 13 UTC to 23 UTC. The atmospheric ther- 
modynamics structure captured by the radiosonde data 
collected at the SGP site at 08 and 23 UTC (Figures lc and 
Id) indicates several features related to the aforementioned 
precipitation activity. A capping inversion of temperature 
was observed under freezing (0°C) level with humid 
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SGP central facility (Cl) site on 25 April 2011. Balloon-borne 


boundary layer. As a result, shallow warm clouds were most 
likely to be formed under the conditions limiting vertical 
development of clouds. However, cold deep layer in condi- 
tional instability is observed above 0°C level in the free 
atmosphere. Deep convections thus tended to be easily 
formed when once strong convection is triggered through the 
capping inversion. 

3. Results 

3.1. Surface Rainfall Structure: Comparison Between 
the Disdrometer Measurements and Simulation 

[7] Laser optical Parsivel disdrometers [Loffler-Mang and 
Joss, 2000] were installed at 16 locations around the ARM 
SGP central site (36.36°N, 97.29°W) (Figure SI). The max- 
imum distance between these disdrometers was 12.26 km. 
The disdrometers had a sampling of the maximum length of 
the falling particle in a horizontal plane and sorted particles 
into a matrix of 32 size intervals ranging 0 to 26.78 mm and 
32 fall velocity intervals ranging from 0.05 to 20.8 m s _1 ; the 
contents in the first and second size intervals (D < 0.26 mm) 
were cut-off because of large noise for small detection. The 


2 of 7 


L24805 


IGUCHI ET AL.: MC3E RAINFALL WRF-SBM ANALYSIS 


L24805 


data was recorded per 10 seconds and integrated to 1 minute. 
The DSD from Parsivel is detemiined from raw 32 x 32 size 
versus fall velocity matrix and is expressed as a function of 
drop diameter, D (mm): 

N ( Dj ) ~ At ' v(Dj)AreaADj ’ ® 

where j is the index of the diameter bin, At is observation 
period in seconds, Area is the cross-sectional area, which is 
nominally 54 cm 2 since length and width of the laser beam are 
180 mm and 30 mm, respectively. But it is expressed as a 
function of D (length x [width — (DJ 2)]) to account partially 
observed drops that are at the edge of the laser beam. AD, is 
the bin width, C ) is the number of drops at a given size bin and 
measured fall velocity, v(D / ). 

[8] Integral rainfall parameters, i.e., total number con- 
centration, N, (drops m -3 ), bulk raindrop effective radius, 
r e (/im), rainfall rate, R (mm hr -1 ), radar reflectivity, 
Z (mm 6 m -3 ), are calculated by integrating DSD spectra as 
follows: 


1 v — >\ 

N, = — V 

A f 


Cj 


At y—f v(Dj)Area ’ 
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32 

Z = J2 D j 6 N(D j )AOj. (5) 

7=1 

Figures 2a and 2c show scatter diagrams of N,(D > 0.77 mm)- 
R and r e -R correlations derived from the 15 Parsivel mea- 
surements (the data of one unit was corrupted) for 24 hours 
on 25 April 2011. Note that N t is calculated using parts of 
DSDs with D > 0.77 mm in consideration of possible 
underestimation of number concentration at the small size 
range in Parsivel measurements [Tokay et al., 2001], The 
plots are color-coded for the period from 06 to 12 UTC in 
blue and from 03 to 06 and from 12 to 24 UTC in red; this 
grouping is based on the difference in approximate time 
of passages of the deep convective system and shallow 
boundary layer clouds. The groups of the blue- and red- 
color dots can be recognized as two distinct modes charac- 
terized by large r e and R and small r e and R, respectively 
(Figure 2c). The threshold value of r e is approximately 
400 /im, and R ranges from 0.01 to 10 mm hr -1 and 
from 0.1 to 100 nun hr -1 in each group. On the other 
hand, Figure 2a shows that N, group of the blue-color 
marks branches off with slightly smaller N, from the red- 
color group at approximately R = 1 mm hr - , while both 
modes has a monotonous increase in R at log-scale with 
a similar increase in N, up to 1000 in 3 . This behavior is 


related to difference in r e . In the red-color group, small 
and almost constant r e suppresses an increase in R heavier 
than 1 0 mm hr - 1 . In contrast, R increases up to 1 00 mm hr - 1 
with an increase in r e in the blue-color group. The existence 
of the two distinct modes for correlation between particle size 
and rainfall rate can be justified by the 2-dimensional video 
disdrometer (2DVD) measurements for the same rainfall 
event (Figure S2b). 

[9] In order to obtain the corresponding precipitation 
DSDs by the disdrometric measurement and subsequently 
integral rainfall parameters from the WRF-SBM simulation, 
so-called Parsivel simulator was designed: Momentary 43 
binned PSDs of the WRF-SBM simulation were recorded 
per hour for all model grid points. Then, surface 32-binned 
precipitation DSD (equation (1)) was calculated using the 
43-binned PSDs and terminal fall velocities for mass bins in 
the SBM. The binned size ranges of the 32 bins were iden- 
tical to those used in the Parsivel disdrometer program. The 
contents in the first and second bins ( D < 0.26 mm) were 
cut-off to be aligned with the actual disdrometer sampling. 
Then, the corresponding integral rainfall parameters were 
calculated following equations (l)-(5). This approach allows 
a reasonable straightforward comparison of rainfall para- 
meters between the simulation and measurement, though 
uncertainty due to the incompatibility of sampling timing 
and location remains. The integral parameters were sampled 
over a domain of 20 km square centered at the ARM SGP 
site. The application of longer distance than the maximum 
distance between the distrometers is to fill up the difference 
in number of sampling between the measurements and 
simulation with different sampling intervals. 

[ 10 ] Figures 2b and 2d show scatter diagrams of N t (D > 
0.77 mm)-R and r e -R correlations derived from the WRF- 
SBM simulation. The plots are color-coded for the period 
from 06 to 12 UTC in blue and from 03 to 05 and from 13 to 
24 UTC in red (the reason of the different notations between 
the measurements and simulation is because the momentary 
output of the simulation was recorded per hour). The simu- 
lation reproduced overall two distinct characteristics of N r R 
and r e -R that correspond reasonably well with those in the 
observed plots: Branch off of the blue-color marks group 
from the red-color group is simulated at R > 1 mm hr -1 
(Figure 2b). Two distinct modes characterized by large 
r e and R (blue) and small r e and R (red) are simulated 
(Figure 2d). However, the simulated plots pattern are more 
scattered than that of the measurements. In particular, dif- 
ference in r e between the measurements and simulation is 
highlighted at R < 1 mm hr - 1 ; the simulated r e is scattered to 
larger value in both color groups. The overprediction of r e at 
the small R range is caused by a group of DSD with small 
N t less than 0.5 m , seen only in the simulation plot 
(Figure S3). The marks with overprediction of r e can be 
almost removed by sampling only DSDs with N, larger 
than 0.5 m -3 (Figure S4). Although the distrometers may 
capture DSDs with the tiny N„ they are disregarded as 
noise [Tokay et al., 2005]. The set of large r e , small N, 
and R is simulated at grids probably where rainfall begins 
to fall and only large raindrops reach to the surface in 
advance. Spontaneous breakup of large drops is not con- 
sidered in the present SBM, so that huge raindrops after 
melting of huge ice particles like hailstones may artifi- 
cially reach to the ground without breakup and be sampled 
in the simulation. 
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a) Parsivel Measurement 
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Figure 2. Scatter diagrams (a, b) between N, and R, (c, d) between r e and R derived from the disdrometric observational 
data and from the WRF-SBM simulation result, and (e, f) between Z and R, derived from the disdrometric observational data 
and from the WRF-SBM simulation result. Plots are color-coded by the sampling period, 06-12 UTC (blue) and 03-06, 12- 
24 UTC (red) for the measurements and 03-05, 13-24 UTC (red) for the simulation. Lines and equations in Figures 2e and 
2f show the fitted lines approximating a power function and the approximate equations for each color group, respectively. 


[n] The relationship between Z and R is generally 
approximated to be a form of a power law function, Z = AR h . 
The values of the coefficients are recognized as an important 
aspect of the rainfall characteristics, though the values have 
non-negligible dependency on the types of measurement 
instruments and the regression techniques [ Campos and 
Zawadzki, 2000; Tokay et al., 2001]. Figures 2e and 2f 
illustrate the Z-R relationships in the form of scatter dia- 
grams derived from the measurements and simulation, 
respectively. The parameters, A and b, are computed using a 
linear regression method separately for the plot groups of 


red-color and blue-color, shown in the form of Z = AR h in 
the panels. The values of A and b of the simulation are 
generally good agreement with those of the observation, 
except for A for the blue-colored dots. This difference is 
considered to be due to the same source as in the afore- 
mentioned difference in r e at the range of R < 1 mm hr - 1 . 

[ 12 ] The timing of appearances of convective and strati- 
form clouds around the SGP site in the simulation is justified 
by the overall agreement between the observation and sim- 
ulation panels in Figure 2, though the present grouping 
according to the time ranges may be coarse. Shallow 
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a) WRF-SBM surface precipitation re (w/o diameter cutoff) [um] 



b) WRF-SBM cloud top temperature [degree C] 



Figure 3. Snapshots of horizontal distributions of (a) surface precipitation r e (without diameter cutoff) and (b) cloud-top 
temperature of the WRF-SBM simulation at (left) 09 and (right) 21 UTC on 25 April 2011. 


boundary clouds with rainfall, however, seem to be over- 
simulated and included in the 06-12 UTC group as com- 
pared to the measurements, because a non-neglectable part 
of blue-color marks are overlapped by the red-color marks in 
the simulation panels of Figure 2. 

3.2. Relationship Between the Surface Rainfall 
Structures and the Cloud/Weather Conditions 

[13] These distinct rainfall characteristics originated from 
difference in processes between the deep convective and 
shallow boundary layer clouds to form the precipitation. 
Figure 3 illustrates the horizontal distributions of r e (without 
cutoff at particular diameter for surface DSD) and cloud-top 
temperature, derived from the WRF-SBM simulation at 09 
and 21 UTC; the cloud-top temperature is here defined as air 
temperature at the highest model level containing the total 
water content larger than 0.3 g m~ 3 in each grid column. 
This figure clearly demonstrates that r e and the cloud-top 
temperature are inversely correlated. In particular, a bound- 
ary around 400 jim of r e in Figure 3 a is in good agreement 
with a border of 0° Celsius of the cloud-top temperature in 
Figure 3b. These results show that the surface precipitation 
of red- and blue-color groups in Figure 2 mostly originates 
from shallow warm and deep mixed-phase rain processes, 
respectively. 


[14] At 09 UTC (Figure 3, left), couples of deep convec- 
tive system are simulated over the ARM SGP site. The 
coldest cloud-top temperature is approximately — 60°C, and 
the corresponding surface DSD r e is more than 900 jim (up 
to 2000 jim as shown in Figure 2d). Shallow warm clouds 
with the cloud-top temperature over 0°C cover the wide area 
of southwestern segment and move toward the ARM SGP 
site by a southwestern wind. The r e of the shallow system is 
approximately 300 jim. At 21 UTC (right panels), a number 
of patchy shallow clouds remain over area including the 
location of ARM SGP site. The r e is approximately 300- 
400 jim. The simulated cloud top heights over the site at 
09 and 21 UTC are in reasonably agreement with those 
observed by the KAZR (Figure la). The WRF-SBM simu- 
lation poorly reproduced the homogenous distribution of 
shallow stratiform clouds at 21 UTC, which is speculated 
by the KAZR profile around the time. An accurate predic- 
tion of cloud fraction is a difficult problem even in a high- 
resolutional simulation [e.g., Khairoutdinov and Randall , 
2003], and this problem needs to be investigated more in 
future works. 

[15] In the case on 25 April, 0°C temperature level were 
located at a height of approximately 3.5 km, so that warm 
clouds below this level were shallow and raindrops reached 
to the ground from clouds immediately. Terminal fall 
velocity with droplet radius of 400 p is assumed to be 
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approximately 2.6 m s -1 at 1000 hPa level in the model 
microphysics [Khain and Sednev, 1995]. Against this fall 
velocity, maximum updraft velocity in the simulated shallow 
warn clouds is roughly 2 m s _I . Thus, particles that grow 
larger up to the radius of approximately 400 fi m cannot be 
held in the clouds and fall to the ground. The small updraft 
velocity in the shallow clouds is a factor of small surface r e . 
On the other hand, 10 m s -1 or much stronger updraft 
velocities in the deep convective system were simulated, so 
that any size particles in the model can be held in the updraft 
region of the clouds. Thus, much larger particles than those 
from the shallow warm clouds can be generated and finally 
fall to the ground. The surface rainfall rate increases with an 
increase in surface-rainfall r e . 

[16] Figure S5a illustrates the simulated thermodynamic 
fields that cause the fonnation of series of deep convective 
clouds at 09 UTC, which corresponds to those observed by 
the KAZR (Figure la). There are two major flows of warm, 
moist air from the southern part and cold, dry air from the 
northern part on 850 hPa level. Parts of the two major flows 
are converged around 36°N, 98°W. This is approximately 
the center of the deep convective system. A couple of warm 
cores are formed around the coordinates on the 850 hPa 
pressure level. Warm air mass is transported to the level of 
free atmosphere in conditional instability, and induces deep 
convections. On a north-south cross-section of the east-side 
of the system along 97.5°W (Figure S5b), a strong vertical 
wind shear is simulated under 800 hPa level. A low-level 
flow from south rises to a higher level at 35°N by the 
opposite wind, and then air mass diverges at the top of deep 
convection, approximately on 200 hPa level. Deep convec- 
tions are continuously induced along the outflow boundaries 
originating from the convection close behind. The deep 
convective system is conveyed northeastward by the back- 
ground southwestern wind. 

[ 17 ] On the other hand, shallow boundary layer clouds are 
simulated over the southwestern segment at 09 UTC. Noc- 
turnal boundary layer clouds are often observed over this 
area, and the detailed mechanism was investigated in Zhu 
et al. [2001], Turbulence driven by wind shear is an 
important factor to maintain well-mixed boundary layer 
and form clouds, in place of radiative heating in daytime. 
Figure S5c shows that shallow clouds are simulated along 
remarkable wind shear near the ground. Relative humidity 
reaches 100% and condensation occurs under or around 
capping inversion of temperature (Figure S6). Unlike the 
segment shown in Figure S5b, deep convective system are 
not developed, probably because of no convergence pro- 
moting breakup of the capping and subsequently involve- 
ment in free atmosphere with conditional instability. The 
boundary layer structure shown in Figure S6 is similar to 
those observed at ARM SGP site before and after the pas- 
sage of the deep convective system (Figures lc and Id). 
These results imply that the observed soundings at the times 
are under a sufficient condition to form shallow boundary 
layer clouds. 

4. Conclusions 

[is] We have analyzed the Parsivel disdrometer observa- 
tions at the ARM SGP site on 25 April, 2011. The mea- 
surement results show unique distinctive structures of 
correlations between the integral rainfall parameters, which 


can be classified into two groups mostly according to the 
types of clouds causing the precipitation. The group with 
large bulk effective particle size and rainfall rate originates 
from series of deep convective clouds, whereas the inverse is 
tme for the group from shallow boundary layer clouds. The 
WRF-SBM weather prediction simulation could success- 
fully reproduce the rainfall microphysical structures. The 
simulation for high rainfall rate is well, whereas the low rate 
rainfall is simulated relatively with defects; note that the 
measurements also have a large uncertainty at low rainfall 
rate, especially R < 0. 1 mm/hr. The discussion on the surface 
rainfall microphysics was expanded to relevance to the 
weather and cloud conditions on the day, by referring to 
the three-dimensional atmospheric field in the simulation. 
The simulation demonstrated that the horizontal distribu- 
tion of the classified precipitation was well inversely corre- 
lated with that of cloud-top temperature; that is, the distinct 
rainfall structures are derived from the difference between 
the deep mixed-phase and shallow warm rain processes in 
clouds aloft. The characteristic structures of cloud and pre- 
cipitation were most likely to be formed under the weather 
condition on the day. Series of deep convective clouds were 
locally and rapidly developed through meso-scale feedback 
process. Shallow boundary layer clouds were easily formed 
under the atmospheric sounding with humid boundary layer 
under temperature inversion. 

[ 19 ] This paper has presented a possibility of the high- 
resolution NWP simulation coupled with one of the up- 
to-date microphysical schemes, even though this is a single 
case study for one-day precipitation event. This rainfall event 
is set as an important case in the MC3E field campaign dur- 
ing April and May 2011. The three-dimensional output 
including hydrometeor PSDs of the WRF-SBM is planned to 
be provided to the MC3E and GPM science communities to 
support their research, e.g., for testing/developing the GPM 
pre-launch precipitation retrieval algorithm. 
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